Perturbational Treatment of the Gravitational Potential Effect on Binary Black Hole 

Evolution 
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Binary black hole (BBH) systems are usually located in the gravitational potential well formed 
by a massive black hole (BH), which is mostly located in the center of a galaxy. In most existing 
studies, the BBH systems are treated as isolated systems, while the effect of the background is 
ignored. The validity of the approximation is based on the belief that the background gravitational 
field from other sources is extremely weak compared with the strong gravitational field produced by 
the BBH itself during the evolution, and can be neglected in gravitational wave detection. However, 
it is still interesting to check how valid this approximation is. In this work, instead of simulating the 
three-BH problem with a fully relativistic treatment, we use a perturbational scheme to investigate 
the effect of the background gravitational potential on the evolution of a BBH, especially on the 
waveform of its gravitational radiation. Four scenarios are considered including the head-on collision 
and the inspiral-to-merger process of a BBH which is either freefalling towards or circularly orbiting 
around a third large BH. The head-on collision and the circular inspiral are two limits of all possible 
configurations. The existence of the background gravitational potential changes the arrival time 
of the gravitational wavefront of a BBH, prolongs the wavelength, and increases the gravitational 
radiation energy. And most interestingly, the background gravitational potential induces the higher- 
order modes of the gravitational wave of a BBH. These interesting phenomena can be explained by 
the gravitational redshift effect and the change of eccentricity of a BBH's orbit from the background 
gravitational potential. Without further studies, these phenomena could introduce complications or 
even mislead people in the identification of the source of gravitational wave and in distinguishing 
the signatures of an isolated BBH from a BBH in a background gravitational potential. 

PACS numbers: 04.25.Dm, 04.30.Db, 95.30.Sf, 97.60.Lf 



I. INTRODUCTION 



In numerical relativity, Einstein's equation can be 
solved numerically without any approximation or sym- 
metry assumption, with the aid of a supercomputer. Be- 
sides the massive computational cost involved, the sta- 
bility issue is also a nontrivial problem in the numerical 
calculations for Einstein's equation. Breakthroughs in 
2005 and 2006 [l|43| shifted the development of numeri- 
cal relativity to a higher pace. Now many numerical rel- 
ativity groups around the world are capable of evolving 
black hole (BH) systems with high accuracy. Moreover, 
a variety of interesting topics, such as long term gravita- 
tional waves J4|-[l4|. gravitational radiation induced BH 
recoil [15l - l3g |. the estimation of the final BH's mass and 
the spin of binary black hole systems pjll - [3a | , and so on 
|37h42| . have been extensively studied over the past few 
years. Numerical relativity has now become an efficient 
tool in the research of general relativity and astrophysics. 
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Gravitational wave detection is an important task for 
the further development of gravitation and general rel- 
ativity. Currently, t he g round-based laser interferom- 
eters such as LIGO Hf, VIRGO H, GEO600 [H, 
and TAMA [46[ are already up and running, while the 
space-based laser interferometer LISA [47[ is still under 
construction. Progress in detection has made the re- 
quirement for the theoretical prediction of gravitational 
wave signal more and more urgent. Numerical relativ- 
ity plays a key role in both aiding gravitational wave 
data analysis [48| and helping to build template banks 
of the phenomenological waveform [49[ for the most im- 
portant gravitational source — binary black hole (BBH) 
coalescence. 

In most theoretical and numerical calculations, BBHs 
are treated as isolated systems, although stellar-mass 
BBHs usually are located in a galaxy with a super- 
massive BH in the center. The validity of the approx- 
imation is based on the idea that a stellar-mass BBH 
is usually far away from the super-massive BH whose 
gravitational potential is negligible compared with the 
strength of the gravitational field produced by the BBH 
itself. Besides the possible super-massive-BH environ- 
ment, it is possible for stellar-mass BBHs to be affected 



by other stellar BHs and/or stars, especially in globular 
clusters. Although the background potentials in these 
scenarios are usually much weaker than the gravitational 
field of the BBH, it still remains as an interesting ques- 
tion as to how valid the approximation is, and to what 
extent the approximation holds. There have been several 
studies of the three-BH problem with numerical relativ- 
ity from various standpoints i37, 38, 50 ;53| . The results 
all indicate that the dynamics and behavior displayed by 
three BHs are qualitatively different from those of an iso- 
lated BBH. For example, in [50|, numerically generated 
BBH initial data sets and post-Newtonian techniques are 
utilized to show that the presence of a third BH has non- 
negligible relativistic effects on the location of a BBH's 
innermost stable circular orbit (ISCO), an increase in 
merging time, and an amplification of the gravitational 
radiation emitted from the BBH. As we all know, infor- 
mation about a BBH such as the initial/final spin, the 
mass ratio and the orientation can be extracted from a 
detected gravitational wave by comparing the wave train 
with the waveform template in the data bank calculated 
with theoretical models. Without proper consideration 
of the environmental effect on the BBH evolution and 
thus on the distortion of the gravitational wave, the ex- 
tracted information might be incorrect and possibly lead 
to misinterpretation and/or misunderstanding. It is also 
possible that the nonlinearity of general relativity might 
give rise to a qualitatively difference due to a tiny differ- 
ence, just like in other nonlinear systems |54| . It is also 
interesting to seek any possible nonlinear effects in the 
three-BH problem which might be helpful in distinguish- 
ing among different scenarios |55| . 

In this work, instead of the heavily numerical calcula- 
tion of the evolution of a three-BH problem with a fully 
relativistic treatment, we use a perturbational scheme, 
along with numerical simulations, to investigate the ef- 
fect of the background gravitational potential on the evo- 
lution of a BBH, especially on the waveform of the grav- 
itational radiation. As with our previously preliminary 
study [56| , the scenarios we considered include the head- 
on collision and the inspiral-to-merger process of a BBH, 
considered as the two limits in all possible configurations, 
in the cases where the BBH is freefalling towards or cir- 
cularly orbiting around the third large BH. Our results 
show that the existence of the background gravitational 
potential changes the arrival time of the gravitational 
wavefront of the BBH, prolongs the wavelength, and in- 
creases the gravitational radiation energy. And most 
interestingly, the background gravitational potential in- 
duces the higher-order modes of the gravitational wave of 
a BBH. These interesting phenomena can be explained by 
the gravitational redshift effect and the change of eccen- 
tricity of the BBH's orbit from the gravitational poten- 
tial. These phenomena could introduce complications or 
even mislead people in the identification of the source of 
gravitational wave without further studies to distinguish 
the signatures of a BBH with a background gravitational 
potential from the ones of an isolated BBH. 
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FIG. 1: A schematic illustration of a BBH system in the 
gravitational potential of a third massive BH. 



The remainder of the paper is organized as follows: In 
Sec. ini we give a description of the code configuration, 
the gauge condition, the initial data, and especially the 
perturbational method for the outer boundary condition. 
Then we consider the four special cases for the evolution 
of a BBH under the influence of a third large BH, fol- 
lowed by a numerical analysis of these results one by one 
in Sec. IIIII We summarize and discuss the viability of 
the perturbational method and the implications of our 
findings in Sec. lIVI Throughout the paper, the geometric 
units with G = c ~ 1 are used. 



II. NUMERICAL TECHNIQUES 

As a preliminary study of the background effect, we 
investigate the dynamics and the waveforms of a BBH 
system in a potential background which comes from a 
third, distant and massive BH. For this three-body sys- 
tem, we do not simulate its three-body dynamics with a 
fully relativistic treatment as in |53| . Instead, a method 
of perturbation is adopted for the problem since the grav- 
itational field from the third BH to the BBH is weaker 
than the field between the two small BHs in the BBH. 
It is obvious that the effect of the third massive BH is 
stronger if it has larger mass 771,3 or the shorter distance 
R to the binary system, or both. Therefore, we would use 
the ratio of m^/R to denote the strength of the gravita- 
tional background resulting from the third massive BH. 
Figure [1] illustrates the configuration considered in this 
work. 



A. Code configuration and Gauge condition 

Our AMSS-NCKU code for solving Einstein's field 
equations based on the BSSN formalism is updated from 
our previous work [57j . For the time integration, we up- 
dated from the iterative second-order Crank-Nicholson 
method to the fourth-order Runge-Kutta method. We 
also developed our own infrastructure for the implemen- 
tation of the fixed mesh refinement and its parallelization 
in addition to using the GrACE package J5g |. Our new 
infrastructure also provides a "moving box" style mesh 
refinement. The box will move to follow the position of 



the BHs, which arc determined by Q 



part of the extrinsic curvature of each hole, 
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(1) 



We follow closely (59|, |60[ in the development of the new 
infrastructure to deal with the interface between two 
mesh levels. Since three time levels are used for the in- 
terpolation in time, we adopt the scheme described in 
[6l| to prepare the initial data. With GrACE, the time 
step for every mesh level is set to be that the Courant 
factor times the finest spatial resolution. However, with 
our new infrastructure, the time step is set by following 
the recipe in [5^, |60|. Typically, the Courant factor is 
set between 0.25 and 0.5. The results obtained with our 
infrastructure and with GrACE's fixed mesh refinement 
are consistent with each other. 

The choice of gauge condition is essential for a long- 
term stability in the numerical evolution to avoid the 
encounter of singularity and large spatial coordinate 
stretch. The "l-|-log" condition for the lapse [62| and the 
F-driver condition for the shift [GJI , which are called "the 
moving puncture gauge" together, have been successfully 
applied to the long-term stable BBH simulations based 
on the BSSN formulation [2|, y] . In [57[ , we have reviewed 
and studied the parameters and the additional advection 
terms for this type of gauge condition used in the com- 
munity. Therefore, in this work, we adopt 



dto. — XiP-'dja = — 2aA', 



(2) 
(3) 



dtB' -Xsp^djB' = dtr~Xil3^djr-vB\ (4) 

where (/(a), 77, Ai, A2, A3, A4) = (a, 2,1, 1,1, 1). The ini- 
tial gauge conditions are chosen duS a = tp^^, /3* =0, and 
B^ = 0. It has been confirmed in [531 that this configu- 
ration allows stable evolutions. 



B. Multi-puncture initial data with spectral 
method 



In this section, we review the puncture scheme and 
then describe how the multi-puncture BH initial data 
is constructed by the multi-domain spectral method ex- 
tended from the LORENE library [64l.l65|. Here we em- 
phasize that our initial data for the three-BH problem is 
from the numerical calculation solving the full Einstein 
constraint equations. 

With the assumption of conformal flatness and maxi- 
mal slicing in the conformal decomposition of the 3 + 1 
formalism of general relativity (see [66| for instance) , the 
constraint equations are greatly simplified and also de- 
coupled. And thus the momentum constraint allows the 
Bowen-York solution [63] for the conformally trace-free 
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(5) 



where n is the spatial unit vector pointing away from 
the puncture, and P and S correspond respectively to 
the linear momentum and the intrinsic angular momen- 
tum of each hole. A^^ can be linearly superposed for 
multi-hole spacctime with A"-^ = J2a=i ^a' since the mo- 
mentum constraint equation is linear in this case. The 
conformal factor then can be solved from the Hamilto- 
nian constraint as 



D^-iP 



1 



-J' 



' A- 4*J 



(6) 



To deal with the physical singularity of the BH, one sep- 
arates out the singular part 



Ips 



1 + 



\z 



(7) 



from the conformal factor -0, where TOq is the mass pa- 
rameter for each puncture and Tq is the coordinate dis- 
tance from each puncture. Equation ^ is the exact 
multi-hole solution at rest satisfying the flat Laplacian 
equation [gg • The desired regular part u = ip — ips sat- 
isfies the elliptic equation 
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(8) 



The existence and uniqueness of the solution has been 
discussed in [69| . 

Solving elliptic equations is computationally expen- 
sive. The spectral method is emphatically suitable to be 
applied to it for its high precision and fast convergence, 
provided the solution is a smooth function. The spectral 
method has been successfully applied to study the rela- 
tivistic star model ^^. And the LORENE library has 
been developed to provide a multi-shell domain frame- 
work for this type of problem. For the puncture initial 
data, due to the at most C^ property of u at the punc- 
ture as discussed in [6^, the direct application of the 
spectral method would only give a polynomial conver- 
gence. The two-puncture data with spectral method was 
first extensively studied in [7l| with a series of coordinate 
transformations which modify and improve the order of 
smoothness of the solution in the new coordinate. This 
scheme has been widely used in the numerical relativity 
community, though it is not easy to be generalized to a 
multi-puncture scenario. The multi-puncture data was 
studied recently with a multi-grid method in [53 ]. 

Our multi-puncture initial data solver is motivated by 
the earlier work [T^ for the excised BBH initial data. 
In this code we cover on each BH a spherical multi-shell 
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FIG. 2: Local convergence test (at the puncture) for the ini- 
tial data of a single boosted BH (dashed-red) with P = 0.2, 
a spinning BH (dotted green) with S = 0.2, and an equal- 
mass BBH (solid) separated by 6M with P — ±0.2, respec- 
tively. The vertical axis indicates the relative difference be- 
tween the solution u at Nr = A'^ and the one at the next 
higher resolution. For the single BH cases, the angular reso- 
lution {Ng,N^) = {Nr,4) which provides sufficient polar an- 
gular resolution for this axisymmetric test. For the BBH case, 
{Ng,N^) = (13, 16). Both of them displays exponential con- 
vergence at lower resolution and polynomial convergence after 
the difference has dropped below certain value which is de- 
pendent on the choice of domain boundaries. 



domain, and split u ~ '^^Ua (the index a runs over the 
number of the BHs) and the puncture equation (|S]) into 



D^Ur, 
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(9) 



Note that only one of the extrinsic curvature tensors is 
spht in the equation. Therefore, the source term in the 
RHS contributes only near each hole. The use of spheri- 
cal polar coordinates is adequate for solving the equation 
near the punctures. Both Aij and ips are known analyti- 
cally and would be set once and for all in the source term. 
The value of Ua is imposed to be zero at the physical outer 
boundary (at infinity) on the outermost, compactified 
shell. And drUa{r = 0)= is also ensured at the punc- 
tures as the inner boundary condition. These equations 
for each hole are then solved iteratively with the Poisson 
solver in the LORENE library until each successive differ- 
ence of Sua is as small as possible, typically 10^^^. In the 
three-BH scenario, the massive hole is fixed and offers a 
fixed gravitational potential background, therefore there 
is only two Poisson equations that need to be solved in 
each iteration. We have also checked the resulting data 
and found that it is very close to the one from the su- 
perposition of the solution of the small binary and the 
solution of the third massive BH. This is expected since 
the massive BH in our study is distant to the binary. 

The convergence tests shown in Fig. [5] for a single 
boosted BH, a spinning BH, and a BBH display a rapid 



convergence to a high precision. It is clear that the con- 
vergence is exponential with low resolution, but turns to 
be polynomial with high resolution, as expected in j73{ . 
The turning point from an exponential convergence to 
a polynomial convergence mainly depends on the choice 
of the domain boundaries, as well as the order of the 
solution's smoothness. However, it is also known that 
the polynomial convergence at higher resolution depends 
only on the singular structure of the solution. 

The momentum parameter for the quasi-circular bi- 
nary is set according to the fitting curve based on the 
helical Killing vector conditions in [7J]. As we note that, 
in the puncture scheme, only the metric of 3-geometry 
was specified, while the initial lapse and shift are cho- 
sen as the moving puncture gauge condition described in 
the beginning of this section. Nevertheless, it is justified 
in our experience that the choice of initial gauge for the 
puncture data has little effect on the resulting physical 
content, at least within our numerical accuracy. 



C. Boundary condition 

In numerical relativity the Sommerfeld radiation 
boundary condition is widely used and also a good ap- 
proximation for no refiection from the boundary at a 
finite distance. For our BBH simulation with a time- 
independent gravitational background, the field variable 
set Q{t,r) = {(^,7ij — riij,K,Aij,T^} receives the gravi- 
tational contribution from the BBH and from the third 
large BH respectively as 



Q{t, r) = QBBH(i, r) + Qzrdir), 



(10) 



where the first term Qbbh in the RHS includes the dy- 
namical outgoing wave, while the second term Qzrd in the 
RHS serves as a fixed gravitational background which are 
trivially zero except the conformal factor of the third BH, 
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In order to update the variable set Q{t, r) at the outer 
boundary, we first extract Qwan(toi f) at the current time 
step to by subtracting from Q{tQ, r) the Q^tdiT) which is 
fixed in the whole evolution. Then, as usual, the Som- 
merfeld boundary condition is applied [75| to obtain the 
numerical data of QBBH(ii,'r) on the outer boundary for 
the next time step ti. And finally we add Q^rdi''') back 
to the new QBBH(ii, J') to obtain Q{ti, r). For the gauge 
variables, we use the usual Sommerfeld boundary condi- 
tion for stability. 

To justify the appropriateness of the modified bound- 
ary condition, we have performed two simulations with 
different outer boundaries and found their difference on 
the trajectory and gravitational waveform of the BBH 
are ignorable, as shown in Fig. [1^] (Further details are 
explained in Sec. IIIIDJ) . Therefore, our boundary condi- 
tion does not introduce unphysical drift to the BBH and 
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FIG. 3: Fourth-order convergence of the waveforms for the 
head-on collisions of an isolated BBH system. The solid line 
represents the difference of *I'4 between the high and medium 
resolutions. The red-dashed line represents the difference be- 
tween the medium and low resolutions with a convergence 
factor C4. 



the effect observed in this work is indeed resulted from 
the third BH. 



III. NUMERICAL RESULTS 

There is an obvious hierarchical structure in the con- 
sidered system, as shown in Fig. [TJ A small BBH system 
and a third BH, which is distant and much more mas- 
sive, form a bigger binary system. From the point view 
of eccentricity e of a binary, there are two limits of the 
possible configurations existing in the universe : a head- 
on collision (e = 1) and (quasi-)circular orbit (e = 0). 
Here we would like to investigate four possible combina- 
tions of the limiting cases in order to offer some clues for 
understanding this kind of astrophysical system. 

The considered hierarchical structure comprises four 
possible cases. They are (A) the head-on-freefall case: 
a BBH is in the process of a head-on collision while its 
center of mass (CoM) freefalls toward a third BH; (B) 
the head-on-orbiting case: a BBH is in the process of a 
head-on collision while its CoM circularly orbits a third 
BH; (C) the inspiral-freefall case: a BBH is in the process 
of a quasi-circular inspiral while its CoM freefalls toward 
a third BH; (D) the inspiral-orbiting case: a BBH is in 
the process of a quasi-circular inspiral while its CoM cir- 
cularly orbits around a third BH. The numerical results 
of the four scenarios will be presented in the following 
subsections. The total mass of the BBH in these cases is 
set to be M = 1. 

We use the similar grid setup as in our previous work 
j57| . There are five levels of grid for the fixed mesh re- 
finement. The physical boundary is put at 64M. The 
domain takes a cubic shape, and the length of the do- 
main in each level shrinks to one half of the length of 
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FIG. 4: The typical waveforms for the head-on collisions of an 
isolated BBH and for the head-on-freefall case in Sec. IIII Al 
The red-dashed (black-solid) line represents the case with 
(without) the third BH. The (£ = 2,m = 0) mode of *4 
detected at r = 40 is plotted. The gravitational potential of 
the third BH causes mainly shifts in time, decreases in ampli- 
tude, and prolongation of the wavelength in the gravitational 
waveform. 



the preceding level. For the finest level, the domain con- 
tains two boxes: they are bounded by (—2 < a; < 2, 
-2 < y < 2, 0.25 < z < 4.25) and (-2 < a; < 2, 
-2 < y < 2, -4.25 < z < -0.25), and are movable 
with the moving-box technique. Before the study on the 
effect of a third BH, the convergence of our code for an 
isolated BBH in a head-on collision is tested. The fourth 
order convergence of \D'4 is shown in Fig. [3] with the con- 
vergence factor 
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(12) 



where hi = 1/32, /i2 = 1/28 and /13 = 1/24 are the 
finest grid widths for the high, median and low resolution, 
respectively. Some minor disagreement for i > 80 was 
expected to come from the outer boundary which could 
give a second-order error. 

In the study of the BBH's evolution under a gravita- 
tional background, we use the Newman-Penrose scalar 
^1^4 instead of the gravitational strain /i"*"/^ to represent 
the gravitational wave. \I'4 is widely used in numerical 
relativity community as a measurement of gravitational 
waveform. So it is convenient to use ^4 and its spheri- 
cally harmonic decompositions, instead of h^'^ , for com- 
parison with previous works. Furthermore, ^4 can be 
linked to h^^'^ straightforwardly as explained in J4ll7a|. 



A. Head-on-freefall case 

We will consider, in this and the next subsections, the 
BBHs in head-on collisions while their CoMs either free- 
fall toward or circularly orbit around the third massive 
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FIG. 5: (a)Initial physical distance between the small BBH, 
(b)proper merging time from the initial physical distance as 
estimated with Eq. (|13p . and (c)proper distance from the cen- 
ter of the small BBH system to the detector, respectively, 
with respect to the strength of the gravitational potential of 
the third BH. (d) Comparison of the proper merging time, 
which are either measured at the detector(squares) or as the 
direct addition(circles) of (b) the (c). The circles and squares 
indicate the numerical data. The slopes of their linear fitting 
are also given. 



BH. The massive BH has the mass m3 and is located 
at (— i?, 0,0). We will vary the mass ma while fix the 
coordinate distance at i? = 1000, and the ratio m^/R will 
represent the strength of different gravitational potential 
due to the third massive BH [771 . 

In Fig. 21 we compare the gravitational waves from the 
head-on collision BBH affected by a third BH and from an 
isolated binary in our earlier result [53| • The BBH in both 
scenarios has the mass parameters mi = m2 = 0.5 and 
is initially located at (0, 0, ±1.1515), which are about the 
ISCO. We can see from this figure that the gravitational 
potential of the third BH causes a time delay, a slight 
decrease of the amplitude of the (, = 2 mode of ^4 and 
a prolongation of the wavelength in the waveform. We 
also observe at the same time the excitation of higher 
order modes with i > 2. The decrease in the amplitude 
mainly comes from the gravitational redshift caused by 
the third BH. Similar to the nonlinear redshift effects of 
an electromagnetic wave in a medium, the potential of 



the third BH broadens the spectrum of the gravitational 
wave instead of simple linear spectral redshift. 

The time delay comes from three factors: (i) the de- 
layed merger; (ii) the prolonged proper distance from the 
source to the detector and (iii) the change in the coor- 
dinate time. It can be understood that the last one is 
simply caused by the change of time gauge due to the 
existence of the third BH, therefore wc only discuss the 
first two factors in detail. 

Factor (i) comes from the effectively larger proper sep- 
aration, i.e., the physical distance d of the BBH, due 
to the curved spacetime background of the third BH. 
Here the physical distance d is defined as the proper 
length of the shortest line connecting two individual ap- 
parent horizons, which is along the 2-axis in the current 
case, and is explicitly calculable at the initial slice with 
d ^ J y^jij dx^dx^ provided these two apparent horizons 
are found. We plot d in first panel of Fig. [S] with re- 
spect to the potential strength of the third BH. Unlike 
the coordinate merging time of the BBH that can be 
obtained from the simulation, it is ambiguous to eval- 
uate the proper merging time because of the singulari- 
ties in the BBH. Due to our pcrturbational treatment of 
the third BH's gravitational effect and the small velocity 
(compared with speed of light) of each BH in the head-on 
BBH, Newtonian mechanics can provide a rough estimate 
of the proper merging time for a given d by virtue of the 
ordinary equation. 
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where m is the mass of individual BH in the BBH, and 
the physical distance between the small binary and the 
third BH is V. In the last term of the above equation, 
V does not appear in the final expression. However, 
the third BH's gravity still affects the physical distance 
d through its effect on the 3-mctric. We integrate the 
above equation to obtain the estimated proper merging 
time for certain d, and show the results in Fig. [SJd. For 
the factor (ii) , the proper distance between the BBH sys- 
tem and the detector is prolonged due to the potential of 
the third BH. In Fig. [5]:, the prolonged proper distance 
from the BBH to the detector was plotted with respect 
to the potential strength of the third BH. The prolonged 
proper distance should be equal to the proper time as 
the gravitational wave propagates in the speed of light. 
The two effects responsible for the proper time delay, one 
in the merger of head-on collision and the other in the 
propagation of the gravitational wave, can account for 
the numerically detected time delay for the gravitational 
wave. To verify this point, we measure the proper time 
of the highest peak of the gravitational wave arriving the 
detector (squares in Fig. [5ji), and compare it with the 
result from the direct addition of the estimated proper 
merging time and the proper propagation time(circles in 
Fig. [SJi) as mentioned. After a linear fit to get the slope 
as the time delay rates with respect to the extra poten- 
tial, it turns out that the measured proper time delay 
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FIG. 6: Left panel: Waveforms of the head-on-freefall case in Sec. lIII Al for different potential strengths of the third BH. The 
amplitude decreases as the third BH's gravitational potential increases. Right panel: duration for three wavelengths starting 
from the wave peak at around f = 45 to the peak at around t = 90. The wavelength increases as the third BH's gravitational 
potential increases. 
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FIG. 7: Accumulated radiated energy with respect to time for 
different gravitational potential strengths. A stronger gravi- 
tational potential results in more energy radiation. 



rate, 61.47 is very close to the sum of the delay rate of 
the estimated proper merging time and the proper prop- 
agation time, as shown in the Fig. [SJi [T^ . 

As mentioned earlier, the redshift effect due to the 
gravitational potential of the third BH vifill decrease the 
amplitude of ^4's ^ = 2 mode as well as prolonging 
the wavelength of the gravitational waveform. They are 
clearly shown in Fig. [BJ where the stronger gravitational 
potential causes smaller amplitude as in the left panel, 
and the longer durations of three cycles (from the first 
to the third peak in the left panel), i.e., the longer wave- 
length as shown in the right panel. 

We also investigate the gravitational energy radiation 
for this case and show the result in Fig. [71 It shows that 
a stronger gravitational potential results in more energy 
radiation. This can be understood as follows: when the 



small BBH system is located in the gravitational poten- 
tial of a third BH, the potential energy will be put into 
the binary system due to the interaction, and some part 
of this energy is carried out by the gravitational wave. 
Then the radiation energy during the head-on collision is 
enhanced by the potential of the third BH. Campanelli 
et al. [501 have studied the ISCO problem for triple BH 
systems. They concluded that the third BH will enhance 
the gravitational radiation of BBH, which is consistent 
with our numerical result. In our perturbational treat- 
ment, the curved background is fixed and thus there is 
no back-reaction between the BBH and the third BH. In 
a more realistic case, the back-reaction will speed up the 
free-falling process of the BBH toward the third BH. 

Besides distorting the 1 = 2 modes, the existence of the 
third BH also induces higher-order modes, especially the 
£ ~ 3 modes, in the BBH's gravitational radiation. The 
left panel of Fig. [S] shows the waveform of the {£ = 3, 
TO = 1) mode of the gravitational radiation induced in 
the current case with respect to the different strength 
of the gravitational potential. The right panel of Fig. [8| 
shows the amplitudes of the (^ = 3, to = 1) mode of 
the gravitational wave, which demonstrates a nonlinear 
growth with respect to the strength of the third BH's 
gravitational potential. In principle, it is possible for the 
nonlinear growth of the £ = 3 mode to be used as a 
signature to distinguish a threc-BH system from a BBH 
system. However, the phenomenon could also complicate 
the identification of the source of the gravitational wave. 
Even in the Newtonian framework, certain three-body 
system with the same quadrupole wave forms ixiay have 
quite distin guis hable higher-order modes of the gravita- 
tional wave [55| . Without a careful understanding on the 
waveform pattern, the gravitational wave emitted from a 
thrce-BH system, like the cases described in this work, 
could be misidcntificd as one from an unequal-mass BBH 
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FIG. 8: Left panel: Waveforms of the {£ = 3, m — 1) mode of ^4 for different background potential strengths. Right panel: 
The largest amplitude of *I'4,3i showing a nonlinear growth pattern. 
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FIG. 9: Waveforms of Re(r»I'4,2o) for four scenarios mentioned 
in Sec. IIIIBI These results show that the relativistic Doppler 
effect on the behavior of the gravitational waveform is much 
less important than the gravitational redshift effect in the 
head-on-orbiting case. 



system in which the higher-order modes of gravitational 
wave also exist. 



B. Head-on-orbiting case 

In this subsection, we study the head-on collision of a 
BBH while the BBH system moves along a circular orbit 
around a third large BH. The orbiting velocity of the cen- 
ter of mass (CoM) of the BBH is initially set according 
to the Newtonian gravity, i.e., w^ « m^/R. This treat- 
ment is valid since the third massive BH is distant to the 
BBH and the orbit of the BBH is set to be circular. In 
fact, the massive BH are so distant such that by the end 
of merger of the BBH, the CoM displacement of BBH is 



only a little fraction of one big orbit, which makes the 
eccentricity of the big orbit not relevant to our study. 
Given the linear momentum parameters and the position 
parameters of the BHs, we solve the puncture initial data 
as described in Sec. IIIBl For the third large BH, we fix 
the coordinate distance at i? = 1000 as in the previous 
subsection, and change the gravitational potential m^/ R 
by varying the mass parameter 7713. 

There are two effects for a BBH orbiting in the gravi- 
tational potential well of a third large BH: the relativistic 
Doppler effect and the gravitational redshift. In order to 
identify the importance of these two effects on the time 
delay in the gravitational waveform in the current case, 
we compare the waveforms obtained from the following 
four scenarios: (i) the head-on collision of an isolated 
BBH; (ii) the head-on collision of a BBH freefalling to- 
wards a third large BH; (iii) the head-on collision of a 
boosted isolated BBH; (iv) the head-on collision of a BBH 
orbiting around a third large BH. Here the mass param- 
eter of the third BH is set to be m^, = 14.4, the same in 
the cases of (ii) and (iv) , and the boosting velocity in (iii) 
is set to be the same as the orbiting velocity in (iv) , the 
magnitude of which is 0.12, according to v^ « rn^/R and 
R = 1000. The results in Fig.lHlshow that the waveforms 
of cases (i) and (iii) are close to each other. So are the 
waveforms of cases (ii) and (iv). This indicates that the 
gravitational redshift of the third large BH should dom- 
inate the effect on the time delay of the gravitational 
wave in the merging process of a BBH. On the contrary, 
the relativistic Doppler effect has only minor influence 
on the time delay of the waveform. Therefore we will ne- 
glect in the following discussion the relativistic Doppler 
effect due to the small orbiting velocity of the CoM of a 
BBH around a third large BH, and focus mainly on its 
gravitational redshift effect. 

In order to let the BBH move slow enough such that 
it will stay inside the computational domain during the 
merger process, and thus we can keep the necessary ac- 
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FIG. 10: Left panel: Waveforms of the head-on-orbiting case in Sec. HIT 51 for different background potential strengths. They 
are qualitatively similar to Fig. (6] except the smaller decreasing amplitude due to the smaller background potential strength. 
Right panel: coordinate time of the highest peak of ^4,20 arriving the detector for different background potential strengths. 



curacy during the generation of gravitational wave, the 
gravitational potential strength of the third BH is set to 
be smaller than in the previous subsection. The numeri- 
cal results for this case show that all phenomena includ- 
ing the time delay of the waveform, the decrease of the 
wave amplitude, the prolongation of the wavelength, the 
enhancement of the energy radiated, and the excitation 
of the higher-order modes are qualitatively similar to the 
ones in the previous case. Therefore we only show here 
the decrease in amplitude of the gravitational waveform 
and the coordinate time delay of the highest peak of ^4,20 
with respect to different potential strengths in Fig. [101 It 
can be seen in the left panel that the decrease in ampli- 
tude of the waveform is less than that obtained in the 
previous subsection due to the smaller gravitational po- 
tential strength of the third BH. Just as in the previous 
subsection, the orientation of the small BBH barely af- 
fects the result in this case because of the negligible tidal 
effect from the third BH in a head-on collision. 



C. Inspiraling-freefall case 

In this and the next subsections, we will consider the 
effect of a third massive BH on an inspiraling BBH. The 
two BHs of the BBH initially have the the irreducible 
masses mi = 7712 « 0.5 and are located at (0,0, i-j). 
The third BH has mass 7713 = 5 x 10"* and is located at 
(i?,0,0). 

In this subsection, we study the inspiral-to-merger evo- 
lution of a BBH while the CoM of the BBH freely faUs to- 
wards the third large BH. The left panel of Fig. [TT] shows 
the evolution of a typical inspiraling trajectory of the 
BBH, from the initial coordinate separation D = 4.5, af- 
fected by the third large BH located at (i? = 5 x 10^, 0, 0). 
We can see from the plot that the gravitational potential 
of the third large BH leads to an increase in the eccen- 



tricity of the trajectory of the BBH, compared with that 
of an isolated BBH. The right panel of Fig. [TT] shows 
a typical comparison of the gravitational waveforms for 
the BBH inspiraling with and without the influence of 
a third large BH. It can be seen that the gravitational 
potential of the third large BH causes a time advance, 
a prolongation of the wavelength, and a somehow com- 
plicated variation in the amplitude of the gravitational 
waveform. 



As explained in Sec. IIII A[ the time delay for the ex- 
istence of a third large BH in a head-on collision arises 
from three factors: the delayed merger process, the pro- 
longed proper distance, and the change of the coordinate 
time. For an inspiraling process of a BBH, there exists 
another effect, i.e., a change of the eccentricity, as shown 
in the left panel of Fig. [TTJ It is obvious that the change 
of the eccentricity is mainly caused by the tidal effect of 
the third large BH. As a BBH system emits more gravi- 
tational radiation in an elliptical orbit than in a circular 
orbit [73, 113 , the increase of eccentricity is expected to 
expedite the merger process of the BBH. Therefore, the 
time shift in the inspiraling process is the result of compe- 
tition between the increase in eccentricity and the others. 
Here we look at this in more detail. 

From the right panel of Fig. [TT] we can see that the 
waveform is affected by both the time delay and time 
advance phenomena: the spurious radiation part of the 
waveform is delayed, while the burst part is advanced. 
The beginning of the waveform, near t = 50, is a result 
of the spurious radiation hidden in the initial data. Thus, 
this part is independent of the merging process of a BBH. 
Here we can see a time delay because of the prolonged 
proper distance from the BBH to the detector, at r = 40, 
due to the existence of the third BH. In the main part 
of the merging process, starting at the coordinate time 
t w 100, we can clearly see the time advance effect. This 
is due to the increase in eccentricity under the influence 



10 





3 












m3/R = 0.01 




2 


/^^^^^"^ 




1 


a (n^\\ 


N 





y\ v->^ 1 1 




-1 


\^^^// 




-2 


^^^^"^^^SS^^ 




-3 





0< 



-3-2-10123 
Y 



0.08 I > ^ 

^■^^ ' W/0BH3 
0.04 -m3/R = 0.01 




80 100 120 140 160 180 200 
Time (M) 



FIG. 11: Left panel: Typical inspiral trajectories for an isolated inspiraling BBH and for the inspiraling-freefall case in Sec. lIIICl 
The solid curve indicates the trajectory of an isolated BBH and shows a quasi-circular inspiral shape. The (red) dashed curve 
indicates the trajectory of the BBH in the gravitational potential well of the third BH with 7713 = 5 x 10^ which is located 
at (_R = 5 X 10^,0,0). There is a change in the eccentricity of the (red) dashed curved when compared with the solid line. 
Right panel: Real components of r'il 4,^20 for an isolated inspiraling BBH and for the inspiraling-freefall case. The red dashed 
(black solid) line indicates the case with (without) the third BH. The detector is located at r = 40. It can be seen that the 
gravitational potential of the third BH causes a tiny time delay in the spurious part and the effect of a time advance in the 
major part of the waveform. 
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FIG. 12: The Re[r^4,2o] for the inspiraling-freefall case for 
different potential strengths of the third BH. The results show 
the effect of advance in time and the prolongation of the wave- 
length in the burst part due to the third BH's gravitational 
potential. The effect of advance in time is obvious and much 
larger than in the spurious radiation (refer to Fig. 1111) . The 
amplitude first increases and then decreases in the wave-train. 
And the variation becomes bigger as the gravitational poten- 
tial of the third BH becomes larger. This detail is explained 
in the context. 



of the gravitational potential of the third large BH. Ob- 
viously the effect of the eccentricity increase overtakes 
other time delay factors in the competition. 

Besides the effect of time advance, there is another 
feature shown in the waveform. If we take the high- 
est peak in the wave-train of ^4 roughly as the time of 



merger [4| , this divides the waveform in the gravitational 
wave-train into two kinds of different behaviors: Before 
the merger, the wave amplitude of ^4 becomes larger 
as the gravitational potential of the third large BH be- 
comes larger. This is because the increase in eccentricity 
of the BBH's trajectory is the dominant effect in speed- 
ing up the inspiral-to-merger process compared to the 
isolated case. This results in stronger gravitational radi- 
ation. After the merger, the wave amplitude of ^4, turns 
out to become smaller as the gravitational potential of 
the third BH becomes larger. The decrease in the am- 
plitude of the gravitational waveform mainly arises from 
the stronger rcdshift effect due to the existence of the 
third BH. These phenomena are illustrated in the right 
panels of Fig. [TT] and in Fig. [121 

In fact, the inspiral- freefall case has been studied in 
[5(j | as an initial data problem. Based on the result of 
its ISCO initial data, it was conjectured in the work that 
the third BH will (1) increase the terminal amplitude of 
the inspiral gravitational waveform; (2) increase the du- 
ration of the pre-plunge phase; and (3) redshift the fre- 
quency of the gravitational wave during the pre-plunge 
phase. Points (1) and (3) are consistent with our cur- 
rent results, that the radiation power is amplified by the 
third BH before the merging process, and that there is a 
gravitational redshift effect. However, point (2) seems to 
conflict with our time-advance result. Therefore, we need 
to look deeper into this point. After a thorough study, we 
discover that as the initial separation of BBH increases, 
the phase change of the gravitational wave (due to the 
effect of the third large BH) transits from a time-delay 
to a time-advance, as shown in Fig . [T3l Note that the 
ISCO separation case studied in [5(| corresponds to the 
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FIG. 13: Time difference of the inspiral-to-merger process 
between an isolated BBH and the inspiral-freefall case in 
Sec. unci Here the merging time is the coordinate time cor- 
responding to the maximal amphtude of the ^4,20 waveform. 
The time difference ranges from being negative, i.e., time- 
delay, to being positive, i.e., time-advance, as the initial sep- 
aration increases. In this case the turning point is around 
r «4.2. 



leftest data point in this figure, which is time-delay. One 
possible explanation for the result is that, since the ini- 
tial separation is small, the gravitational redshift effect 
caused by the existence of a third large BH is dominant 
in the evolution of the BBH and the merging time is de- 
layed. On the other hand, the eccentricity effect by a 
third large BH outweighs the gravitational redshift effect 
when the initial separation of the BBH becomes large 
enough, and thus the merging time is advanced. In other 
words, the time difference in the inspiral-to-merger pro- 
cess of a BBH under the influence of a third large BH is 
dependent on the initial separation. 

In Fig. [131 We do not expect the time difference to con- 
verge when the separation range is still very short, com- 
pared with a fixed R (the distance to the third BH). Since 
the tidal effect on the evolution of the BBH due to the 
existence of the third large BH is accumulative during the 
BBH's inspiral, we could expect that the time-difference 
will keep increasing as the initial separation of the BBH 



becomes larger [81|. It is known that the gravitational 



waveform from the merger of a BBH with the ISCO as 
the initial separation should agree with the waveforms 
just before the merger from different initial separations. 
However, this understanding is based on the assumption 
that the BBH is isolated. It is not the case when there ex- 
ists an external gravitational potential background. The 
final result should come from the competition between 
the circularization effect from the gravitational radiation 
and the tidal effect from the external gravitational poten- 
tial background. It seems not exist related investigations 
on this issue in the literature. Therefore, we turn to 
check our numerical result. In Fig. [T31 it shows that the 
gravitational waveforms from different initial separations 
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FIG. 14: Time-shifted waveforms of r>I'4,20 in the inspiral-to- 
merger process of a BBH, with different initial separations, 
under the influence of a third BH in the inspiral-freefall case. 
T — Tshift = is the shifted merging time defined at the max- 
imal amplitude, and these waveforms do not match with each 
other despite of the shifting. 



between the two BHs in a BBH system under an exter- 
nal gravitational potential do not match with one an- 
other even the merging times have been shifted to be the 
same. From our numerical experiments, it shows that the 
gravitational circularization effect could not relax totally 
the accumulated eccentricity from the tidal effect caused 
by the external gravitational potential background. A 
detailed investigation is needed to have a better under- 
standing on this phenomenon in the future. 

Compared with the head-on collision of a BBH as dis- 
cussed in Sec. IIII Al the higher-order mode effect due to 
the existence of a third BH in the inspiraling process 
is more evident. The left panel of Fig. [15] shows the 
amplitude of Re(r^4 33) produced during the inspiral- 
to-merger process of a BBH with respect to different 
potential strengths of a third BH. The higher-order 
mode, Re(r^4_33), shows a nonlinear dependence on the 
strength of the gravitational potential, as shown in the 
right panel of Fig. [151 This result is consistent with those 
detailed in the previous two subsections: This nonlinear 
phenomenon could be a signature distinguishing a three- 
BH system from a BBH system, as well as complicating 
the identification of the source of a gravitational wave. 



D. Inspiraling-orbiting case 

In this subsection, we investigate the inspiral-to- 
mcrgcr process of a BBH orbiting around a third large 
BH. Initially, the CoM of the BBH moves around the 
third large BH in a circular orbit. Therefore, the velocity 
of the CoM of the BBH can be approximated by Newto- 
nian mechanics. Similar to the previous subsection, the 
two small BHs with the irreducible mass m\ ~ 7712 ~ 0.5 
are located at (0,0, ±-j) initially. We put the third BH 
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FIG. 15: Left panel: Comparison of Re(r>l/4,33) in the inspiraling-freefall case with respect to different potential strengths of a 
third large BH. The existence of the third BH induces the higher-order {£ = 3) mode and quite a large amount of the radiation 
is emitted away in this mode. Right panel: Highest peaks of Re(r<I'4,33) with respect to different potential strengths of the 
third BH. The nonlinear behavior of the higher-order {£ — 3) mode with respect to the gravitational potential strength of the 
third BH is shown. 
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FIG. 16: Left panel: Typical trajectories for an isolated inspiraling BBH and the inspiraling-orbiting case in Sec. IIIIDI The 
solid quasi-circular inspiral curve is for an isolated BBH. The red dashed curve indicates the trajectory under the gravitational 
potential of the third BH with 7713 = 5 x 10* at (R — 2.5 x 10^,0,0). The blue long-dashed line represents the same case 
but with a doubled outer boundary. Right panel: The corresponding waveforms of r^4,2o. The red-dashed (black-solid) line 
indicates the case with (without) the third BH. The blue long-dashed line represents the case same as the one for the red dashed 
curve but with a doubled outer boundary. It can be seen that, besides the decrease in the wave amplitude, the gravitational 
potential of the third BH causes a tiny time delay in the spurious part and a time advance in the major part of the waveform. 



at (i?, 0,0). In order to get a circular orbit for the BBH 
around the third large BH, we use the quantity y^m^/R 
to evaluate the velocity of the CoM of the BBH with re- 
spect to the third BH. This velocity is then vectorially 
added to the quasi-circular orbiting velocity of each BH of 
the BBH. We have ever tried measuring the eccentricity 
of the small BBH affected by the background to check 
if any extra eccentricity is introduced by the vectorial 
addition. It turns out that there is no observable extra 
eccentricity introduced by this addition during the whole 
inspiral stage. For the third large BH, we fix its mass 
parameter m^ = 5 x 10^ as in the previous subsection 



while varying its coordinate distance R. Given the mass 
parameters, the linear momenta and the positions of the 
three BHs, the puncture initial data was constructed as 
described in Sec. IIIBI 

The left panel of Fig. [16] shows a typical trajectory of 
the BBH with the initial coordinate separation D = 4.5 
for the current case, where the third large BH is located 
at (i? = 2.5 X 10^^,0,0). It can be understood from the 
plot that the CoM of the BBH moves roughly along a 
Newtonian circular orbit around the third large BH. The 
trajectory is in fact the combination of the orbit of the 
CoM of the BBH and the BBH's inspiraling around the 
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FIG. 17: The Re(r>I'4,2o) in the inspiraling-orbiting case with 
respect to different potential strengths of the third BH. Due 
to the existence of the third BH, the time advance effect is 
clear and different from the spurious radiation part, and the 
one for the inspiraling-freefall case (compare Fig. I12|) . The 
amplitude decreases as the gravitational potential of the third 
BH becomes stronger. 



CoM, and it is quite deformed compared virith the trajec- 
tory for an isolated BBH. 

The gravitational waveform of the BBH in the current 
case is presented in the right panel of Fig. [161 Compared 
with the right panel of Fig. [T^ we find that the effect 
of the third large BH on the waveform is stronger in the 
current case than in the inspiraling-freefall cases. In the 
plot, the waveform is delayed in the spurious part but 
advanced in the major part, caused by the gravitational 
potential of the third large BH. This is similar to the 
previous inspiraling-freefall case. In Fig. 1161 we also plot 
the trajectory and ^4 for the same scenario but with 
a doubled outer boundary. The perfect coincidence of 
them indicates that the boundary condition described in 
Sec. Ill Cl is reliable. 

In Fig. [T71 the time advance of Re(r\l/4.2o) due to the 
third large BH is clear and different from its spurious 
radiation part. This phenomenon is similar to the one 
in the inspiraling-freefall case (compared with Fig. [T^ . 
Stronger gravitational potential of the third BH results 
in the smaller amplitude of the waveform. 

The shift of the waveform due to the gravitational 
background of the third BH are quantitatively shown 
in Fig. [THl in which it indicates the transition from the 
time-delay to time-advance as the initial separation of 
BBH increases. Similar to the previous case, the possible 
explanation for the change in the time difference of the 
inspiral-to-merger process of a BBH under the influence 
of a third large BH compared with the one of an isolated 
BBH, arises from the competition between the eccentric- 
ity effect and the gravitational redshift effect. Compared 
with the previous (inspiraling-freefall) case in which the 
(.£ = m = 3) mode of ^4 is dominant among the higher 
modes, the (£ = 3, m = 0) modes of ^4 is dominant in the 
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FIG. 18: Time difference of the inspiral-to-merger process 
between an isolated BBH and the inspiraling-orbiting case in 
Sec. IIIIDI with respect to the initial separation. Here the 
merging time is the coordinate time corresponding to the 
maximal amplitude of the ^4,20 waveform. The time dif- 
ference ranges from being negative, i.e., time-delay, to being 
positive, i.e., time advance, as the initial separation increases. 
In this case, the turning point is around r ~ 3.5. 



current case, as shown in the left panel of Fig. [191 The 
orbiting of the BBH along the x-z plane accounts for this 
effect. In addition, the amplitude of the (£ = 3,to = 0) 
mode of Re(^4) increases nonlinearly. 



IV. SUMMARY AND DISCUSSION 

Motivated by the fact that most BBH systems are lo- 
cated in the gravitational potential of a super-massive 
BH hosted in the center of a galaxy, we investigate the 
effect of the potential on the dynamics of a BBH sys- 
tem. Instead of the heavily numerical calculation of the 
evolution of the three-BH problem with a fully relativis- 
tic treatment, we use a perturbational scheme to investi- 
gate the effect of the gravitational potential from a third 
large BH on the evolution of a BBH, especially on the 
waveform of gravitational radiation. In our perturbation 
method, we ignore the back-reaction of the BBH system 
on the third large BH, regarding the third BH as a back- 
ground for the BBH system. 

The scenarios we consider include the head-on collision 
and the inspiral-to-merger process of a BBH in the cases 
of the BBH system freefalling towards, or circularly or- 
biting around the third large BH, which are considered 
as the two limits in all possible configurations. The ef- 
fect of the gravitational potential from a third BH on a 
BBH system in our study includes: (1) the gravitational 
redshift effect, including the prolongation of proper dis- 
tance, the prolongation of the wavelength, the decrease 
of the waveform amplitude; (2) the increase in the eccen- 
tricity which is from the tidal effect of the third large BH 
and expedites the merger of a BBH system; (3) delaying 
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FIG. 19: Left panel: Comparison of Re(r^4,3o) in the inspiraling-orbiting case for different potential strengths of a third BH 
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to different potential strengths of the third BH. The nonlinear behavior of the I = Z mode and the behavior of the curve are 
different from the one in the inspiraling-freefall case. 



or advancing the merging process of the BBH, which de- 
pends on the competition between the gravitational red- 
shift effect and the eccentricity effect; (4) inducing the 
higher-order modes in the gravitational waveform. The 
orientation of the BBH system hardly affects the above 
results. It is interesting that (4) supports the conjecture 
proposed in [5a] in the full GR regime. This might pro- 
vide valuable information for gravitational wave astron- 
omy to distinguish whether the detected BBH system is 
isolated or located in the potential of a super-massive 
BH. 

In this work, for the ^ < 3 modes, all the quantities are 
linearly dependent on the gravitational potential m^/R, 
which validates the usage of the perturbation method. In 
most realistic cases, the gravitational potential is much 
weaker than the cases considered in this work. Therefore 
the perturbation method is expected to be applicable to 
most realistic cases. 

The third BH in our case introduces a non- 
axisymmetric background feature, which makes the issue 
of gravitational wave extraction even more complicated. 
However, considering the detector located at the weak 
field region and finitely separated from both the merging 
binary and the massive BH, we adopt the usual approach 
to extract the information of gravitational radiation via 
the Newman-Penrose Scalar 'I'4. A possible extension of 
the work would be a detailed study on the more rigorous 
way to extract the radiation under a background and the 
dependence of the modification on the relative position 



of the binary, detector and the background source. 

We have considered the evolution of an equal-mass 
BBH under the infiuence of a third large BH. The study 
can be generalized to an unequal-mass BBH system. The 
major effects from the third large BH on a BBH are 
shown qualitatively in this study. The results indicate 
that these phenomena could introduce complications or 
even be misleading in the identification of the source of 
gravitational wave without further study to distinguish 
the signatures of a BBH in a background gravitational 
potential from the ones of an isolated BBH. 
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